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ABSTRACT 

We analyze all X-ray timing data on IE 1207.4—5209 in supernova remnant PKS 1209—51/52 
gathered in 2000-2005, and find a highly stable rotation with P = 424.130751(4) ms and P = 
(9.6±9.4) x 10~ 17 s s _1 . This refutes previous claims of large timing irregularities in these data. In the 
dipole spin-down formalism, the 2a upper limit on P implies an energy loss rate E < 1.5 XlO 32 ergs s _1 , 
surface magnetic field strength B p < 3.5 x 10 11 G, and characteristic age r c = P/2P > 24 Myr. This 
t c exceeds the remnant age by 3 orders of magnitude, requiring that the pulsar was born spinning 
at its present period. The X-ray luminosity of IE 1207.4-5209, L hoi « 2 x 10 33 (d/2 kpc) 2 ergs s" 1 , 
exceeds its E, implying that Lboi derives from residual cooling, and perhaps partly from accretion of 
supernova debris. The upper limit on B p is small enough to favor the electron cyclotron model for at 
least one of the prominent absorption lines in its soft X-ray spectrum. This is the second demonstrable 
case of a pulsar born spinning slowly and with a weak B-field, after PSR J1852+0040 in Kesteven 79. 
Subject headings: ISM: individual (PKS 1209-51/52) — pulsars: individual (IE 1207.4-5209, PSR 
J1852+0040) — stars: neutron — supernova remnants 



1. INTRODUCTION 

The neutron star IE 1207.4-5209 in the center of su- 
pernova remnant PKS 1209—51/52 is the first discovered 
(Helfand & Becker 1984) and most intensively studied of 
the so-called Central Compact Objects (CCOs). These 
seemingly isolated NSs are defined by their steady flux, 
predominantly thermal X-ray emission, lack of optical or 
radio counterparts, and absence of a surrounding pul- 
sar wind nebula (see Pavlov et al. 2004, for a review). 
IE 1207.4—5209 acquired special importance when it be- 
came the first CCO in which pulsations were detected 
(Zavlin et al. 2000; Pavlov et al. 2002). It was distin- 
guished again as the first isolated NS to display strong 
absorption lines in its X-ray spectrum (Sanwal et al. 
2002; Mereghetti et al. 2002; Bignami et al. 2003). 

More recently, accumulated X-ray observations of 
IE 1207.4—5209 were presented as showing large- 
amplitude changes of both sign in its spin period (Zavlin 
et al. 2004) that were unlike any other pulsar and difficult 
to explain. Consequently, the surface dipole magnetic 
field, which is a key parameter in all proposed mecha- 
nisms for the X-ray absorption lines, could not be es- 
timated independently from the spin-down rate, which 
was indeterminate. In this Letter, we present a defini- 
tive study of the spin history of IE 1207.4-5209 that 
corrects previous errors in the data and their analysis. 
We provide reliable spin parameters and discuss their im- 
plications for the interpretation of the X-ray spectrum of 
IE 1207.4-5209, and for the origin of the class of CCOs 
more generally. 

2. ARCHIVAL X-RAY OBSERVATIONS (2000-2005) 

We reanalyzed all timing data on IE 1207.4—5209 from 
the archives of the Newton X-Ray Multi-Mirror Mission 
(XMM-Newton) and Chandra observatories. They span 
2000 January to 2005 July. A log of these observations 
is given in Table 1. 



All 11 XMM-Newton observations of IE 1207.4-5209 
used the pn detector of the European Photon Imaging 
Camera (EPIC-pn) in "small window" (SW) mode to 
achieve 5.7 ms time resolution. Several EPIC-pn data 
sets had photon timing errors uncorrected in their orig- 
inal processing (Kirsch et al. 2004). We reprocessed 
all EPIC data using the emchain and epchain scripts 
under Science Analysis System (SAS) version xmm- 
sas_20060628_1801-7.0.0, which produces correct photon 
time assignments. The observations were affected by 
background to varying degree. To maximize the signal- 
to-noise ratio in each, we adjusted the source extraction 
aperture individually. For this soft source, an energy cut 
of 0.5 — 2.5 keV was found to maximize pulsed power. 

We also examined data simultaneously available from 
the EPIC MOS camera, operated in "full frame" mode. 
Although not useful for timing purposes (2.7 s readout), 
the location of the source at the center of the on-axis 
MOS CCD allows a better background measurement to 
test for flux variability, an important indicator of accre- 
tion, than the EPIC-pn SW mode. The seven observa- 
tions of 2005 exhibit root-mean-square source flux vari- 
ability of less than 1% over the 40 day span. Comparing 
these count rates to the earlier XMM-Newton observa- 
tion of 2001 December, implies a marginally significant 
flux decrease of 5% ± 3% during the 5 year interval. 

Four Chandra observations suitable for timing mea- 
surements of IE 1207.4—5209 are available. They used 
the Advanced Camera for Imaging and Spectroscopy 
(ACIS) in continuous-clocking (CC) mode to provide 
time resolution of 2.85 ms. Two of the four observa- 
tions were taken with the LETG transmission grating in 
place, the zeroth order image being used for timing. This 
study uses data processed by the latest pipeline software 
(revision v7.6.8.1), with the exception of 2000 January 
6, which is processed with revision v6.5.1. Reduction 
and analysis used the standard software package CIAO 
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Log of X-ray Timing Observations and Summary of Results 
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Set 3 combined; 
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Chandra 


ACIS-S/CC b 
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2003 Jun 10 
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Set 5 combined: 


3393.6 


53543.515 


424.1307512(40) 


411.9 



a Period derived from a Z\ test. Uncertainty in last digits is in parenthesis, which is la computed by the Monte Carlo method described 
in Gotthelf ct al. (1999). 
b These Chandra observations used the low-energy transmission grating (LETG). 



(v3.4) and CALDB (v3.3). The photon arrival times in 
CC mode are adjusted in the standard processing to ac- 
count for the known position of the pulsar, spacecraft 
dither, and SIM offset. These needed corrections were a 
potential cause of timing errors in earlier work, but are 
now accurately performed. 

3. TIMING ANALYSIS 

For each observation in Table 1 we transformed the 
photon arrival times from IE 1207.4—5209 to Barycen- 
tric Dynamical Time (TDB) using the coordinates given 
in Table 2, and identified the pulsed signal using a stan- 
dard FFT. To obtain a most precise value of the period 
in each observation we then generated a periodogram 
using a Z\ test (Buccheri et al. 1983) around the FFT 
value. The la uncertainty in P was determined by the 
Monte-Carlo method described in Gotthelf et al. (1999). 
In contrast to the large period changes claimed by Za- 
vlin et al. (2004), a linear ephemcris is an excellent fit to 
our derived periods (see Fig. 1), with xl = 0.65 for 12 
degrees-of-freedom, and no significant detection of a pe- 
riod derivative. The average period throughout the data 
span is P — 424.130801(57) ms, with a formal 2<r upper 

limit of P < 3.9 x 10~ 15 s s" 1 . 

In order to increase the precision, we refitted closely 
spaced observations coherently wherever possible. The 
results for these combined data are listed by set in Ta- 
ble 1, and the methods are described here. The 2005 
June- July set of seven observations spanned 40 days 
specifically to obtain a phase-coherent timing solution. 
Woods et al. (2006) were not able to eliminate large phase 
residuals between these observations or find a unique so- 
lution. We determined that the original data processing 
had timing errors. 

Using the reprocessed data, we iteratively measured 
the period and phase of adjacent observations of the 2005 
data by the Z\ method. Starting with the the longest ob- 
servation of 2005 July 11, we extrapolated the resulting 



period to the flanking observations, verifying that the 
predicted phase and its uncertainty matched to < 0.1 
cycles the actual phase derived from the adjacent obser- 
vations. After completing this procedure for all the 2005 
observations a coherent fit was obtained for the entire 
set. Figure 2 shows the phase residuals of the individual 
observations from the best fit, which demonstrates the 
validity of the solution. The best-fit period has Z\ = 412 
and agrees with the period found above from the inco- 
herent analysis of all observations, while the next highest 
peak in the power spectrum has Zf = 275 and is clearly 
an alias. To test for a P, we then performed a Z\ search 
on a two-dimensional grid of P and P. The extra param- 
eter did not increase the peak Z\ significantly, meaning 
no detection of P. 




500 1000 1500 2000 

TIME (DAYS SINCE MJD 51500) 

Fig. 1. — Period residuals after fitting a linear solution to individ- 
ual observations (top panel) and grouped data sets (bottom panel) 
from Table 1. The error bar for the first individual observation is 
used, and continues off-scale, in the bottom panel. 



The two XMM-Newton observations of 2002 August 4 
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TABLE 2 
Spin Parameters of IE 1207.4- 



5209 



Parameter 


Value 


Right ascension, R.A. (J2000) a 


12 h 10 m 00 s .91 


Declination, Decl. (J2000) a 


-52°26'28".4 


Epoch (MJD) 


53562 


Spin period, P (s) 


0.424130751(4) 


Period derivative, P 


(9.6 ±9.4) x 10~ 17 


Valid range of dates (MJD) 


51549-53582 


Surface dipole magnetic field, B p (G) b 


< 3.5 x 10 11 


Spin-down luminosity, E (ergs s~ 1 ) b 


< 1.5 x 10 32 


Characteristic age, t c (Myr) b 


> 24 



a Mcasurcd from Chandra ACIS-I ObsID 3913, in agreement with 
Wang ct al. (2007). 
b Quantity derived from 2<x upper limit on P. 



and 6 were easily joined, resulting in the period listed 
in Table 1. Now knowing the precise and consistent 
values of P in 2002 and 2005, we were able to make 
a phase-connected combination of the 2001 December 
XMM-Newton observation and the 2002 January Chan- 
dra one, which arc 13 days apart, by finding an exact 
period match for the correct peak from among nearby 
aliases. Finally, we combined the set of two Chan- 
dra observations spanning 2003 June 10-19, which again 
yielded a consistent period at the highest peak in the Z\ 
periodogram. After making these coherent combinations, 
it was not possible to achieve a further phase-connected 
solution over a longer time span, as the intervening cy- 
cle counts could not be determined uniquely. Therefore, 
we made a final linear least-squares fit to the five points 
listed in Table 1, yielding the ephemeris presented in Ta- 
ble 2. This piecewise coherent measurement, yielding a 
2(7 upper limit of P < 2.8 x 10~ 16 , is more than an or- 
der of magnitude more precise than the fully incoherent 
analysis (see Fig. 1). 




-10 10 

TIME (DAYS SINCE MJD 53563.148) 

Fig. 2. — Pulse phase residuals for the 2005 XMM-Newton obser- 
vations of IE 1207.4—5209 after fitting a coherent timing solution 
using a constant period model, yielding P = 424.1307512(40) ms. 



4. INTERPRETATION 

Contrary to previous claims, the timing behavior of 
IE 1207.4—5209 does not require glitches, a binary com- 
panion, and perhaps not even accretion of fallback mate- 
rial, although the latter may still be needed to contribute 
to its X-ray spectrum and luminosity. The absence of de- 
tectable spin variations is due to a weak dipole magnetic 
field. In the dipole spin-down formalism, the 2a upper 
limit on P implies, for an isolated pulsar, an energy loss 



rate E = -inn = 4n 2 IP/P 3 < 1.5 x 10 32 ergs s~ 

19, 



surface magnetic field strength B p = 3.2 x 10 9 V PP < 
3.5 x 10 11 G, and characteristic age r c = P/2P > 24 Myr. 
In its spin properties, IE 1207.4—5209 is nearly a twin 
of another CCO, PSR J1852+0040 (Gotthelf et al. 2005; 
Halpern et al. 2007). The next section closely follows the 
discussion in Halpern ct al. (2007) , which anticipated the 
present result. 

4.1. Cooling and/or Accreting 

The X-ray luminosity of IE 1207.4-5209 is of thermal 
origin, with L bo \ w 2x 10 33 (d/2 kpc) 2 ergs s -1 (De Luca 
et al. 2004). This is much larger than the upper limit 
on its spin-down power, E, and argues that it is mostly 
residual cooling. However, fits to the spectrum require 
two blackbody components; the hotter one, of temper- 
ature A;Tbb = 0.32 keV, has an area of only 0.87 km 2 
(Dc Luca et al. 2004), which may indicate heating by 
accretion onto the polar cap. The canonical area of the 
open-field-line polar cap is R pc = 2w 2 R 3 /Pc « 0.27 km 2 . 
This is ~ 30% of the fitted blackbody component, but 
accretion may cover a wider area. 

The characteristic age r c > 24 Myr, compared to the 
remnant age, estimated as 7 kyr with an uncertainty of 
a factor of 3 (Roger et al. 1988), requires that pulsar was 
born spinning at its current period. A recent popula- 
tion analysis of radio pulsars favors a wide distribution of 
birth periods (Faucher-Giguere & Kaspi 2006), in which 
424 ms would be typical. Furthermore, as magnetic 
field is generated by a turbulent dynamo whose strength 
depends on the rotation rate of the proto-neutron star 
(Thompson & Duncan 1993), it is natural that pulsars 
born spinning slowly would have the weaker S-fields; the 
model of Bonanno et al. (2006) supports this. 

There are no young radio pulsars with B p < 10 11 G. 
Since IE 1207.4—5209 is not necessarily beyond the ra- 
dio pulsar death line, either empirical (Faucher-Giguere 
& Kaspi 2006) or theoretical (Chen & Ruderman 1993), 
there may be another reason it is radio quiet. It is 
possible that low-level accretion of SN debris prevents 
CCOs from becoming radio pulsars for thousands or even 
millions of years. Accretion from a fallback disk (Al- 
par 2001; Shi & Xu 2003; Ek§i et al. 2005; Liu et al. 
2006) was one of the theories considered by Zavlin et al. 
(2004) to explain the now defunct timing irregularities 
of IE 1207.4-5209. But accretion may still be needed to 
account for its radio-quiet and X-ray-hot properties. 

If the magnetic field is weak enough that an accretion 
disk can penetrate the light cylinder, the hotter portion 
of the NS surface in IE 1207.4-5209 can be powered 
by accretion of only m s=s 10 13 g s _1 . It is required that 
B p <5x 10 11 G for IE 1207.4-5209 to be able to accrete 
in the propeller regime. But in this limit, the pulsar 
would tend to spin down at a rate that is excluded by 
observations, 



1.2 x 10" 



14 ,8/7^3/7 



M29 



13 



M 



-2/7 



£ 45 



P 



0.424 s 



where the magnetic moment /i = B p R 3 /2 w 10 29 S Pi n 
G cm 3 . Also, M would have to be greater than 10 13 
g s _1 , as most of the accreting matter is expelled from 
the magnetospheric radius rather than accreted. 
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However, if B p < 2 x 10 9 G, then IE 1207.4-5209 may 
accrete as a "slow rotator," and spin up at a small rate, 



-2.1 x 10- 17 i4 7 m% 7 



M \ 



3/7 



0.424 s 



In this regime, secular spin-up, and torque noise, which 
may be of the same magnitude, are below the sensitivity 
of the existing measurements. 

While flickering is also an indicator of accretion, we do 
not have strong evidence of variability of IE 1207.4—5209 
(< 1% on month timescales). Also, upper limits on opti- 
cal/IR emission from IE 1207.4—5209 are comparable to 
that expected from a geometrically thin, optically thick 
disk accreting at the rate required to account for its X- 
ray luminosity (Zavlin et al. 2004; Wang et al. 2007). 
Therefore, it may be necessary to invoke a radiatively 
inefficient flow in order to consider accretion. 

4.2. X-ray Absorption Lines 

Broad absorption lines in the soft X-ray spectrum of 
IE 1207.4-5209 are centered at 0.7 keV and 1.4 keV 
(Sanwal et al. 2002; Mcrcghctti ct al. 2002), and possibly 
at 2.1 keV and 2.8 keV (Bignami et al. 2003; De Luca et 
al. 2004), although the reality of the two higher-energy 
features has been disputed (Mori et al. 2005). Proposed 
absorption mechanisms include electron cyclotron in a 
weak (8 x 10 10 G) magnetic field (Bignami ct al. 2003; 
De Luca et al. 2004), atomic features from singly ionized 
helium in a strong (2 x 10 14 G) field (Sanwal et al. 2002; 
Pavlov & Bezchastnov 2005), and iron (Mcrcghctti et al. 
2002), or oxygen/neon in a normal (10 12 G) field (Hailcy 
& Mori 2002; Mori & Hailey 2006). 

Our upper limit, B p < 3.5 x 10 11 G, favors the electron 
cyclotron model, for at least one of the lines, over all oth- 
ers that require stronger fields. The cyclotron prediction, 
8 x 10 10 G, assumes that 0.7 keV is the fundamental en- 
ergy E c = 1.16(B/10 11 G)/(l + z), where z is the gravita- 
tional redshift. Another solution postulates hydrogenic 
oxygen for the 0.7 keV kine, while the 1.4 keV line is 
the cyclotron fundamental (Hailey & Mori 2002; Mori & 
Hailcy 2006). As these authors pointed out, abundant 
oxygen may be accreted from supernova debris. One 
caveat, however, is that the magnetic field strength at 
the NS surface can be larger in places than the global 
dipole that determines the spin-down rate. 

4.3. Are CCOs a Class? 

The half dozen radio-quiet CCOs are similar in their 
X-ray luminosities, high temperatures, and absence of 
pulsar wind nebulae. Therefore, they may comprise a 
fairly uniform class defined by a weak magnetic field, 
which in turn results from slow natal rotation. If accret- 
ing, a slow initial spin is still unavoidable, since even at 



Eddington-limited accretion rates, the spin-up and spin- 
down time scales in §3.1 and in Halpern et al. (2007) are 
much longer than the ages of the remnants. While we 
do not have definite evidence of accretion in any CCO, 
small Bp and large P both make it possible for a pulsar 
to accrete at low rates from a SN debris disk; the large 
B p and rapid spin of young radio pulsars prevents such 
a disk from penetrating the light cylinder. In order to 
test the general applicability of these results to the class 
of CCOs, more sensitives searches for their pulsations 
are required. However, prior null results on all of them 
suggest that their pulsed amplitudes are very small. 



5. CONCLUSIONS 

A comprehensive analysis of all timing data on the X- 
ray pulsar IE 1207.4—5209 has resolved the dilemma of 
its mysterious spin properties by correcting previous er- 
rors in data processing and analysis. It is simply a low 
magnetic field NS that has no discernible variation in 
spin over 5 years. If an isolated NS, the upper limit 
on its spin-down power is much less than its bolomct- 
ric X-ray luminosity, which leaves only internal cool- 
ing and/or accretion as possible energy sources. In ci- 
ther case, IE 1207.4—5209 must have been born with a 
weak magnetic field and its long rotation period. We 
speculate that these two parameters are causally related 
and, projecting from the near twins IE 1207.4—5209 and 
PSR J1852+0040, possibly the physical basis of the CCO 
class. 

An additional benefit of this solution is the new 
constraint on proposed absorption-line models for 
IE 1207.4-5209 that depend on the magnetic field 
strength. Our upper limit on the dipole field is close to 
the prediction of the electron cyclotron model for both 
lines, and perhaps oxygen for one of the lines. To actu- 
ally make a significant measurement of B p as small as 
8 x 10 10 G from dipole spin-down would require a fully 
phase-coherent timing solution spanning > 6 yr, assum- 
ing that there is no glitch or other timing noise. Such a 
program would also be sensitive to accretion torques at 
the lowest rates predicted here for spin-up. In case B p 
is as small as 10 9 G and some of the X-ray luminosity of 
IE 1207.4—5209 is due to accretion from a debris disk, 
spin-up could be detected. 



This investigation is based on observations obtained 
with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA Member 
States and NASA. Support for this work was provided by 
NASA through XMM grant NNX06AH95G and Chandra 
Award SAO GO6-7048X. 
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